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Abstract 

The dynamics in QCD with different number of fermion flavors Nf 
is discussed. The emphasis is on the description of the conformal phase 
transition with respect to Nf separating a phase with no confinement and 
chiral symmetry breaking and a phase with confinement and with chiral 
symmetry breaking. 

1 Introduction 

The infrared dynamics of QCD can be viewed only dimly via presently available 
tools. Because of that, there may be surprises, as it has already happened with 
our understanding of nonperturbative dynamics in N—l supersymmetric QCD 
[Q. In particular, those studies showed that the conventional wisdom which has 
accepted that the asymptotic freedom implies confinement is not always true. 

Recently, there has been considerable interest in the existence of a nontrivial 
conformal dynamics in 3+1 dimensional non-supersymmetric vector like gauge 
theories, with a relatively large number of fermion flavors Nf || [| 0, ||, [| 0, || . 
The roots of this problem go back to a work of Banks and Zaks || who were 
first to discuss the consequences of the existence of an infrared-stable fixed point 
a = a* for Nf > N5 in vector-like gauge theories. The value N* depends on 
the gauge group: in the case of SU(3) gauge group, Nf = 8 in the two-loop 
approximation. 

A new insight in this problem ^ || has been, on the one hand, connected 
with using the results of the analysis of the Schwinger-Dyson (SD) equations 
describing chiral symmetry breaking in QCD (for a review, see Refs.@ 0) 
and, on the other hand, with the discovery of the conformal window in N = 1 
supersymmetric QCD [Q. 
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In particular, Appelquist, Terning, and Wijewardhana Q suggested that, in 
the case of the gauge group S\J(N C ), the critical value Nj r ~ AN C separates a 
phase with no confinement and chiral symmetry breaking (Nf > N^ r ) and a 
phase with confinement and with chiral symmetry breaking (Nf < Nj r ). The 
basic point for this suggestion was the observation that at Nf > Nj r the value 
of the infrared fixed point a* is smaller than a critical value a cr ~ n^~i f ' 
presumably needed to generate the chiral condensate |l0|, pj| . 

The authors of Ref. Q considered only the case when the running coupling 
constant a(/x) is less than the fixed point a*. In this case the dynamics is 
asymptotically free (at short distances) both at Nf < NT and N^ r < Nf < 

Yamawaki and the author [|3| analyzed the dynamics in the whole (a,Nf) 
plane and suggested the (a, iV/)-phase diagram of the SU(iV c ) theory (see Fig. 1 
below) .Q In particular, it was pointed out that one can get an interesting non- 
asymptotically free dynamics when the bare coupling constant a<°) is larger 
than a*, though not very large. 

The dynamics with > a* admits a continuum limit and is interesting 
in itself. Also, its better understanding can be important for establishing the 
conformal window in lattice computer simulations of the SU(iV c ) theory with 
such large values of Nf. In order to illustrate this, let us consider the following 
example. For N c = 3 and Nf = 16, the value of the infrared fixed point a* is 
small: a* ~0.04 (see below). To reach the asymptotically free phase, one needs 
to take the bare coupling less than this value of a*. However, because of 
large finite size effects, the lattice computer simulations of the SU(3) theory with 
such a small would be unreliable. Therefore, in this case, it is necessary to 
consider the dynamics with a(fi) > a*. 

The existence of the phase transition with respect to Nf in QCD raises the 
following question: what are the infrared properties of the QCD (3 function 
for different Nf and how the (3 function structure reflects the existence of this 
phase transition? I will discuss those issues at the end of my talk but first I will 
discuss the dynamics in the conformal window of QCD in detail. In particular, 
I will consider the spectrum of low energy excitations in that dynamics in the 
presence of a bare fermion mass Jl^] . We will see that in this case, unlike the 
familiar QCD with a small Nf (Nf=2 or 3), glueballs are much lighter than 
bound states composed of fermions, if the value of the infrared fixed point is 
not too large. Another characteristic point is a strong (and simple) dependence 
of the masses of all the colorless bound states on the bare fermion mass, even if 
the latter is tiny. 

This talk is based on papers []l2|, [l3| . 



1 This phase diagram is essentially different from the original Banks-Zaks diagram Hi. For 
details, see Sec. VII in Ref. J^] 
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2 Dynamics in the Conformal Window in QCD 
like Theories 



I begin by recalling the basic facts concerning the two-loop function in an 
SU(iV c ) theory. The function is 



= -ba 2 - co? (1) 



with 14 



b = ±-(UN c -2N f ), (2) 
1 A^ 2 — 1 

While these two coefficients are invariant under change of a renormalization 
scheme, the higher-order coefficients are scheme dependent. Actually, there is 
a renormalization scheme in which the two-loop function is (perturbatively) 
exact |n|. We will use such a renormalization scheme. 

If b > (Nf < Nf* = ^y^) and c < 0, the function has a zero, corre- 
sponding to a infrared-stable fixed point, at 

a = a* = - h -. (4) 

When Nf is close to NJ*, the value of a* is small. For example, from Eqs.(||), 
(|), and @, one gets a* ~ 0.04, 0.14, 0.28, and 0.47 for N c =3 and N f =W, 15, 
14, and 13, respectively. 

The value of a* becomes equal to a cr = f a t c l° se to Nf ~ 4iV c . 
And the fixed point disappears at the value Nf — Nf, when the coefficient c 
becomes positive (Nt is Nf ~ 8.05 for N c =3). 

The function (|lj) leads to the following solution for the running coupling: 

5 log (l) =— 1 1 / a(g)(a(n)-a') \ 

\fij a(q) a((i) a* \a(ji)(a(q) - a*) J ' 

We emphasize that this solution is valid both for < a* and a(/t) > a* . 

Let us first consider the case with a(p) < a*. It is convenient to introduce 
the parameter [1] 



A = (i exp 



1 fa* — a(u) \ 1 



ba* \ a(/i) ) ba{p) 
Then, Eqs. (g) and (|) imply that 



(6) 



1 -^(?)4^fAi. 



a(<?) VA/ a* \a* — a(q) 
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Taking q = A, wc find that 



1 + e- 



0.73a* < a (A) < a* 



(8) 



One may think that A plays here the same role as Aqcd in the confinement 
phase. However, as we will see, its physical meaning is somewhat different. 
Eq. (^) implies that 

for q >> A (the usual behavior in asymptotically free theories), and 

a(q) 



1 



(10) 



for <7 << A, governed by the infrared fixed point a*. 

Let us turn to a less familiar case with a(/i) > a*. One still can use Eq.(||). 
Introduce now the parameter A as 



A = fi exp 



a(/i) — a* 



ba(fi) 



(compare with Eq.(^)). Then, Eqs.(^) and ( [Tl| ) imply 



a(q) 



&log — 
A 



1 

a* 



log 



a(q) 



a(q) — a* 



What is the meaning of A? It is a Landau pole at which ct(q)\ ^ 
taking q — A in Eq. (|T^) , one gets 

1 1 , a(A) 
— =- = — log — = . 

a (A) a* a(A)-a* 



(11) 



(12) 



Indeed, 



(13) 



The only solution of this equation is a (A) = oo. 

The presence of the Landau pole implies that the dynamics is not asymptoti- 
cally free. To get a more insight in this dynamics, let us introduce an ultraviolet 
cutoff M with the bare coupling constant = a(q)\ q= M- Now all momenta 
q satisfy q < M. 

Eq.(12) implies that at finite — a(M), the cutoff M is less than A, with 
a(A) — oo. Therefore the Landau pole is unreachable in the theory with cut 
off M and with a' ' < oo. Still one can of course use A (11) for a convenient 
parametrization of the running coupling a(q) (see Eq.fll^)). However, one should 
remember that momenta q satisfy q < M < A. 



Eq.(12) implies that 



a\q) 



26 log ^ 



(14) 
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for a(q) » a*, and 

* l- e - a i(^' (15) 

when a(q) is close to a*, i.e. when — a* << a*. Thus, now a(g) approaches 
the fixed point a* from above (compare with Eq.(^0|)). And, in general, Eq.(ll) 
implies that a(q) monotonically decrease with q, from a (q) = a(°' at q = M to 
cc(g) = a* at q = 0. 

Does a meaningful continuum limit exist in this case? The answer is of 
course "y es "- As it follows from Eq.jl^), when M (and therefore A) goes to 
infinity, and the bare coupling > a* is arbitrary but fixed, a(q) is equal to 
the fixed value, a(g) — a*, for all q < oo. Therefore it is a non-trivial conformal 
field theory. 

So far we considered the solution for a(q) connected with the perturbative 
(and perturbatively exact in the 't Hooft renormalization scheme |l5| ) f3 function 
(|j). However, unlike ultraviolet stable fixed points, defining dynamics at high 
momenta, infrared-stable fixed points (defining dynamics at low momenta) are 
very sensitive to nonperturbative dynamics leading to the generation of particle 
masses. For example, if fermions acquire a dynamical mass, they decouple from 
the infrared dynamics, and therefore the perturbative infrared fixed point (|]) 
will disappear. 

The phase diagram in the (a^°\ iV/)-plane in this theory was suggested in 
Ref.||. It is shown in Fig. 1. The left-hand portion of the curve in this figure 
coincides with the line of the infrared-stable fixed points a*(Nf) in Eq.(||). It 
separates two symmetric phases, Sx and S2, with < a* and > a *, 
respectively. Its lower end is N f = Nf (with Nf ~ AN C if a cr ~ 7^5; f ): 
at < Nf < N"' the infrared fixed point is washed out by generating a 
dynamical fermion mass. 

The horizontal, Nf — Nf r , line describes a phase transition between the sym- 
metric phase Si and the phase with confinement and chiral symmetry breaking. 
As it was suggested in Refs.|^|, ||, based on a similarity of this phase transi- 
tion with that in quenched QED 4 0, [ll], |l|] and in QED 3 0, there is the 
following scaling law for m^ yn : 



mlyn ~ A^exp ( ? == = ] (16) 



C 

'WJ) 

where the constant C is of order one and A cr is a scale at which the running 
coupling is of order a cr . 

It is a continuous phase transition with an essential singularity at Nf = Nf r . 
The characteristic point of this phase transition is that the critical line Nf = 
Nr r separates phases with essentially different spectra of low energy excitations 

H D and the different structure of the equation for the divergence of the 
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dilatation current (i.e. with essentially different realizations of the conformal 
symmetry) Q . It was called the conformal phase transition in Ref . [|| . 

At present it is still unclear whether the phase transition on the line Nf = 
Nj r is indeed a continuous phase transition with an essential singularity or it is 
a first order phase transition [Q |7j. However, anyway, the two properties (the 
abrupt change of the spectrum of excitations and the different structure of the 
equation for the divergence of the dilatation current in those two phases) have 
to take place. 

At last, the right-hand portion of the curve on the diagram occurs because at 
large enough values of the bare coupling, spontaneous chiral symmetry breaking 
takes place for any number Nf of fermion flavors. This portion describes a phase 
transition called a bulk phase transition in the literature, and it is presumably a 
first order phase transition. ^ The vertical line ends above Nf=0 since in pure 
gluodynamics there is apparently no phase transition between weak-coupling 
and strong-coupling phases. 

Up to now we have considered the case of a chiral invariant action. But how 
will the dynamics change if a bare fermion mass term is added in the action? 
This question is in particular relevant for lattice computer simulations: for 
studying a chiral phase transition on a finite lattice, it is necessary to introduce 
a bare fermion mass. We will show that adding even an arbitrary small bare 
fermion mass results in a dramatic changing the dynamics both in the S\ and 
S*2 phases. 

Recall that in the case of confinement SU(A C ) theories, with a small, Nf < 
Nf r , number of fermion flavors, the role of a bare fermion mass is minor 
if m' ' << Aqcd (where Aqcd is a confinement scale). The only relevant 
consequence is that massless Nambu-Goldstone pseudoscalars get a small mass 
(the PC AC dynamics). 

The reason for that is the fact that the scale Aqcd , connected with a scale 
anomaly, describes the breakdown of the conformal symmetry connected both 
with perturbative and nonperturbative dynamics: the running coupling and 
the formation of bound state. Certainly, a small bare mass « Aqcd is 
irrelevant for the dynamics of those bound states. 

Now let us turn to the phase Si and 5*2, with Nf > Nf T . At finite A in S\ 
and A in 52, there is still conformal anomaly: because of the running of the 
effective coupling constant, the conformal symmetry is broken. It is restored 
only if A — > in Si and A > M — > oo in 5*2 . However, the essential difference 
with respect to confinement theories is that both A and A have nothing with 
the dynamics forming bound states: since at Nf > Nj r the effective coupling 

2 The fact that spontaneous chiral symmetry breaking takes place for any number of fermion 
flavors, if a(°) is large enough, is valid at least for lattice theories with Kogut-Susskind 
fermions. Notice however that since the bulk phase transition is a lattice artifact, the form of 
this portion of the curve can depend on the type of fermions used in simulations (for details, 
see Ref. §). 
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is relatively weak, it is impossible to form bound states from massless fermions 
and gluons (recall that the Si and 5*2 phases are chiral invariant). 

Therefore the absence of a mass for fermions and gluons is a key point for 
not creating bound states in those phases. The situation changes dramatically if 
a bare fermion mass is introduced: indeed, even weak gauge, Coulomb-like, in- 
teractions can easily produce bound states composed of massive constituents, as 
it happens, for example, in QED, where electron-positron (positronium) bound 
states are present. 

To be concrete, let us first consider the case when all fermions have the same 
bare mass It leads to a mass function m(q 2 ) = B(q 2 )/A(q 2 ) in the fermion 

propagator G(q) = (qA(q 2 ) — B(q 2 ))~ 1 . The current fermion mass m is given 
by the relation 

m{q 2 )\ q i =m i = m. (17) 

For the clearest exposition, let us consider a particular theory with a finite 
cutoff M and the bare coupling constant — a(q)\ q= M being not far away 
from the fixed point a*. Then, the mass function is changing in the "walking" 
regime |L8| with a(q 2 ) ~ a*. It is 



m{q 2 ) ~ m<°> (jY" (18) 

where the anomalous dimension 7 m ~ 1 — (1 — .2— )V 2 0, 0- Eqs.(|i~7|) and 
@ imply that 

rn^^L)^. (19) 

There are two main consequences of the presence of the bare mass: 

(a) bound states, composed of fermions, occur in the spectrum of the theory. 
The mass of a n-body bound state is ~ nm; 

(b) At momenta q < m, fermions and their bound states decouple. There 
is a pure SU(A r c ) Yang-Mills theory with confinement. Its spectrum contains 
glueballs. 

To estimate glueball masses, notice that at momenta q < m, the running of 
the coupling is defined by the parameter b of the Yang-Mills theory, 

b=^N c . (20) 

Therefore the glueball masses M g i are of order 

Aym ^ mexp(- J^). (21) 



ba 

For N c = 3, we find from Eqs.(§), (1), and @ that exp(-^i r ) is 6 x 10 

If 



2 x 10~ 2 , 10 -1 , and 3 x 10" 1 for N f =16, 15, 14, and 13, respectively. Therefore 
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at Nf =16, 15 and 14, the glueball masses are essentially lighter than the masses 
of the bound states composed of fermions. The situation is similar to that in 
confinement QCD with heavy quarks, m >> Aqcd- However, there is now a 
new important point: in the conformal window, any value of (and therefore 
m) is "heavy": the fermion mass m sets a new scale in the theory, and the 
confinement scale Aym (0) is less, and rather often much less, than this scale 
m. 

This leads to a spectacular "experimental" signature of the conformal win- 
dow in lattice computer simulations: glueball masses rapidly, as (mW) i+~t m , 
decrease with the bare fermion mass for all values of less than cutoff 
M. 

Few comments are in order: 

(1) The phases <Si and S2 have essentially the same long distance dynamics. 
They are distinguished only by their dynamics at short distances: while the 
dynamics of the phase Si is asymptotically free, that of the phase S% is not. In 
particular, when all fermions are massive (with the current mass m), the con- 
tinuum limit M — > 00 of the ^-theory is a non-asymptotically free confinement 
theory. Its spectrum includes colorless bound states composed of fermions and 
gluons. For q < m the running coupling a(q) is the same as in pure SU(iV c ) 
Yang-Mills theory, and for all q > m a(q) is very close to a* ("walking", ac- 
tually, "standing" dynamics). For those values Nf for which a* is small (as 
Nf=16, 15 and 14 at N C =S), glueballs are much lighter than the bound states 
composed of fermions. Notice that, unlike the case with m = 0, there exists an 
S-matrix in this theory. 

(2) In order to get the clearest exposition, we assumed such estimates as 

Nf ~ AN C for Nf and j m = 1 — ^Jl — -^f- for the anomalous dimension 

7 m . While the latter should be reasonable for a* < a cr (and especially for 
a* << a cr ) [|l0| 0, the former is based on the assumption that a cr ~ jf^rjf 
which, though seems reasonable, might be crude for some values of N c . It is 
clear however that the dynamical picture presented in this paper is essentially 
independent of those assumptions. 

(3) So far we have considered the case when all fermions have the same bare 
mass m,(°\ The generalization to the case when different fermions may have 
different bare masses is evident. 

(4) Lattice computer simulations of the SU(3) theory with a relatively large 
number of Nf [[l9[ ^o) indeed indicate on the existence of a symmetric phase. 

However, the value of the critical number Nj r is different in different simu- 
lations: it varies from Nf = 6 @ through Nf = 12 |§. 

I hope that the signature of the conformal window considered in this talk 
can be useful to settle this important issue. 
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3 Pade- Summation Approach to QCD /^-function 
Infrared Properties 

How does the structure of the QCD (3 function reflect the existence of the phase 
transition with respect to Ntl This problem has been addressed in the work 

s 

The impetus for that work was the structure of the (3 function of the SU (N c ) 
SUSY gluodynamics which is known exactly if no matter fields are present plf : 



P( X > = 7~ 



1 



1 - N c x/2 



x=-. (22) 

7T 



This structure shows the following two noticeable features: 

a) Eq. (|22|) implies that, besides the conventional asymptotically free phase 
with x < 2/N c , there exists a strong ultraviolet phase in which the coupling x 
is greater than the (3- function pole at x = 2/N c (p2|. In that phase, the value 
x = oo is an ultraviolet fixed point. As in the case of the infrared fixed point 
considered in the previous section, these two phases share common infrared 
properties. Because of that, the /^-function pole is called an infrared attractor. 
Unlike theories with an infrared fixed point, theories with an infrared attractor 
correspond to the confinement phase. 

In Ref. we addressed whether Pade-summations of the MS QCD (3 

function for a given number of flavors exhibit an infrared fixed point, or alter- 
natively, an infrared attractor. The main results are the following. Below an 
approximant-dependent flavor threshold (6 < Nf < 8), the Pade-summation 
(3 functions incorporating [2|1], [1|2], [2|2], [1|3], and [3|1] approximants whose 
Maclaurin expansions match known higher-than-one-loop contributions to the 
/3-function series always exhibit a positive pole prior to the occurrence of their 
first positive zero, precluding any identification of this first positive zero as an 
infrared fixed point. This result is shown to be true regardless of the magnitude 
of the presently-unknown five-loop /3-function contribution explicitly appearing 
within Pade-summation f3 functions incorporating [2|2], [1|3], and [3|1] approx- 
imants. Like in the case of supersymmetric gluodynamics, the pole in question 
suggests the occurrence of dynamics in which both a strong and an asymptot- 
ically free phase share a common infrared attractor. As Nf increases above an 
approximant-dependent flavor threshold, Pade-summation (3 functions exhibit 
dynamics controlled by an infrared fixed point. This fixed point decreases in 
magnitude with increasing flavor number. 

Thus utilizing Pade-summation QCD (3 functions, we obtain a good degree of 
agreement with infrared properties predicted j^, ^, ^| via the 't Hooft renormal- 
ization scheme Jl5| in which the (3 function is truncated subsequent to two-loop 
order. It is noticeable that the infrared structure of the (3 function we obtained 
in QCD with a small number of flavors is similar to that of the [3 function in 
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SUSY gluodynamics: there are strong arguments in the literature in the support 
of essentially the same mechanism of confinement in those theories. 
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Figure 1: The phase diagram in an SU(iV c ) gauge model. The coupling con- 
stant and S and A denote symmetric and asymmetric phases, 
respectively. 
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